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[i]  This  study  is  in  an  area  south  of  the  Shandong  peninsula,  near  the  region  where  Zhou 
et  al.  (1991)  observed  anomalous  drops  in  acoustical  intensity.  Solitary  wave  generation 
and  propagation  simulations  are  performed  using  the  Lamb  (1994)  nonhydrostatic 
model.  The  model  simulations  show  that,  for  summer  conditions,  the  existing  semi-diurnal 
tidal  flow  over  the  topographic  variations  formed  internal  bores  and  solitary  waves.  For 
the  Shandong  area,  we  analyzed  summer  observations  from  Synthetic  Aperture  Radar 
(SAR)  that  tracked  solitary  wave  trains  from  their  surface  roughness  signatures.  The 
images  contained  seven  events  consisting  of  internal  bores  and  solitary  waves  that  traveled 
in  a  well-defined  direction  for  2.5  days.  The  origin  of  the  trains  appeared  at  a  well-defined 
point  along  a  steep  topographic  drop.  The  SAR  observations  guided  and  tuned  the 
model  simulations,  by  comparing  spectra  of  observed  and  modeled  wavelengths.  The 
tuned  model  yields  wavelengths  within  factors  of  2,  or  less,  of  those  derived  from  SAR 
data.  Wavelength  and  amplitude  dispersion  analysis  showed  two  dispersion  regimes. 

Modeled  phase  speeds  were  at  the  lower  limit  of  phase  speeds  deduced  from  SAR  data, 
from  about  0.8  to  1.0  m/s.  Acoustical  intensity  calculations  in  the  presence  of  solitary 
wave  trains  will  be  undertaken  in  a  subsequent  paper  using  a  parabolic  equation  acoustical 
model  along  the  path  of  solitary  wave  train  propagation. 
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1.  Introduction 

[2]  Recent  interest  in  the  oceanography  of  a  region  of  the 
Yellow  Sea  south  of  the  Shandong  Peninsula  can  be 
attributed  to  the  acoustical  measurements  of  shallow-water 
sound  propagation  published  by  Zhou  et  al.  [1991].  They 
made  acoustical  measurements  in  this  region  over  a  period 
of  years.  Some  of  their  measurements  showed  an  anomalous 
drop  in  acoustical  intensity  of  about  15-20  dB  at  a  range  of 
28  km  from  the  acoustic  source.  These  anomalous  losses 
occurred  over  a  narrow  band  of  acoustic  frequencies.  For 
the  case  they  published,  the  loss  was  centered  around  630  Hz. 
The  anomalous  losses  were  both  seasonally  and  directionally 
dependent,  occurring  only  in  the  summer  and  along  tracks 
perpendicular  to  the  shore.  The  authors  postulated  the  exis¬ 
tence  of  solitary  waves  in  the  thermocline,  and  attributed  the 
anomalous  losses  to  the  interactions  of  the  acoustic  signal 
with  the  solitary  internal  waves.  The  lack  of  oceanographic 
data  required  them  to  use  a  model  to  simulate  the  solitary 
wave  environment;  a  gated  sine  function  was  used  to  repre¬ 
sent  the  solitons.  A  parabolic  equation  acoustic  propagation 
model  was  used  to  perform  the  transmission  loss  calculations. 
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The  simulation  results  from  this  hypothetical  environment 
showed  that  an  anomalous  transmission  loss  could  occur  near 
the  same  frequency  (630  Hz)  as  evident  in  the  acoustic 
measurements.  Computer  simulations  by  other  investigators 
[Chin-Bing  et  al .,  1993;  King  et  al .,  1994]  subsequently 
confirmed  that  a  “resonant-like”  transmission  loss  is  caused 
by  an  acoustical  mode  coupling  in  the  presence  of  solitary 
waves,  together  with  a  corresponding  larger  bottom  attenua¬ 
tion  for  the  coupled  acoustic  modes.  However,  the  lack  of 
supporting  oceanographic  information  left  some  doubt  as  to 
the  soliton  hypothesis.  Schooling  fish  [ Weston ,  1992]  were 
suggested  as  a  possible  source  of  the  anomalous  loss  seen  in 
the  data. 

[3]  In  the  acoustic  calculations  the  existence  of  solitary 
waves  has  so  far  been  only  postulated  for  the  area  south  of 
the  Shandong  Peninsula.  This  paper  addresses  this  issue  by 
considering  solitary  wave  generation  and  propagation  in  the 
region. 

[4]  The  generation  and  propagation  of  solitary  waves 
have  been  modeled  with  nonhydrostatic  models.  Nonhydro¬ 
static  modeling  of  tidal  flow  over  steep  topography  shows 
that  tidal  flow  depresses  the  pycnocline  and  generates  a 
depression  [ Lamb ,  1994;  Brandt  et  al.,  1997;  Wam-Varnas 
et  al .,  2003].  From  the  depression,  two  internal  bores  are 
formed  propagating  in  opposite  directions.  The  offshore 
internal  bores  propagate  and  steepen  on  the  leading  edge 
through  nonlinear  effects.  Then  frequency  and  amplitude 
dispersion  sets  in  and  the  leading  edge  disintegrates  into 
solitary  waves  [Lamb,  1994;  Brandt  et  al .,  1997;  Warn- 
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Figure  la.  Track  of  the  RADARSAT1  measurements,  black  line,  located  south  of  the  Shandong 
peninsula.  The  smaller  black  lines  reflect  the  orientation  of  the  internal  bores  and  solitary  wave  trains 
relative  to  the  direction  of  propagation. 


Vanias  et  a l .,  2003].  The  onshore  propagating  internal  bore 
develops  a  hydraulic  jump  during  semidiurnal  tidal  reversal, 
when  the  flow  becomes  supercritical  [ Lamb ,  1994;  Brand! 
et  ai. ,  1997;  Warn-Varnas  et  ai ,  2003]. 

[5]  The  paper  is  organized  into  the  following  sections.  In 
section  2  the  Shandong  region  and  SAR  observations  of 
solitary  waves  are  described.  In  section  3  the  numerical 
simulations  of  the  generation  and  propagation  of  solitary 
wave  trains  are  presented  with  a  discussion  of  the  results. 
We  Finish  with  conclusions  in  section  4.  Simulations  of  the 
acoustical  field  structure  and  intensity  will  be  published  in  a 
separate  paper. 

2.  Shandong  Region 

2.1.  Topography  of  the  Region 

[6]  The  Yellow  Sea  is  a  shallow  sea  bounded  by  China, 
Korea,  and  the  East  China  Sea  with  a  maximum  depth  of 


around  80  m  in  the  center.  The  topography  varies  gently 
from  the  Chinese  continent  and  more  rapidly  from  the 
Korean  peninsula.  The  present  study  is  located  south  of 
the  Shandong  peninsula  and  will  be  referred  to  as  the 
Shandong  area.  Figure  la.  The  black  line  shows  the  track 
of  a  solitary  wave  train  observed  with  RADARS  ATI 
Synthetic  Aperture  Radar  (SAR).  At  the  beginning  of  the 
arrow,  there  is  a  relatively  steeper  slope,  where  the  first 
internal  bore  is  observed.  The  lines  across  the  arrow  indicate 
the  along  crest  direction  of  the  wave  packets.  The  variable 
angle  of  the  lines  suggests  refraction  along  the  shelf  break. 

[7]  The  topography  of  the  area  is  extracted  from  the 
Digital  Atlas  of  Choi  [1999].  This  is  a  bathymetric  database 
for  the  Yellow  Sea  and  neighboring  Sea  of  the  Korean 
Peninsula.  Its  construction  encompassed  all  available  ba¬ 
thymetry  charts.  A  detailed  view  of  the  three-dimensional 
topography  in  the  area  where  the  satellite  transects  are 
located  is  shown  in  Figure  lb.  A  mountainous  type  of 
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Figure  lb.  Three-dimensional  topographic  distribution 
according  to  the  digital  atlas  of  Choi  [1999].  The  solid 
black  lines  represent  40-m  and  70-m  tracks  from  the  shelf  to 
deeper  areas.  The  70-m  track  is  the  same  as  in  Figure  la. 
Vertical  scale  is  the  location  relative  to  the  shelf  depth  in 
meters.  Horizontal  scales  are  in  kilometers. 


distribution  is  evident.  Along  both  of  the  transects,  there  is  a 
drop  nearshore,  followed  by  a  rise  and  a  second  drop.  The 
maximum  depth  of  the  drops  depends  on  which  transect  is 
used.  The  two  transects  used,  indicated  in  Figure  I,  have 


final  depths  of  about  40  and  70  m.  In  Figure  2  the 
topography  along  one  transect,  along  with  the  model 
representations,  is  shown.  The  topography  is  represented 
as  a  combination  of  hyperbolic  functions  (shown  as  the 
dashed  line  in  Figure  2)  because  a  continuous  representation 
of  the  bottom  is  more  convenient  for  the  simulations. 

2.2.  SAR  Data 

[s]  Synthetic  Aperture  Radar  detects  the  variations  in  sea 
surface  roughness  caused  by  internal  wave  currents.  The 
convergence  and  divergence  of  the  internal  wave  currents 
causes  surface  roughness  variations  that  translates  into 
bright  and  dark  patterns  of  pixel  intensity  in  the  SAR 
images. 

[9]  Summer  SAR  observations  were  obtained,  for  the 
Shandong  area,  from  RADARSAT1  SAR  with  a  500  by 
500  km  width  observation  area.  The  observations  were 
acquired  on  8  August  1998  and  processed  at  the  Alaska 
SAR  Facility.  The  pixel  spacing  is  100  m.  Figure  la  shows 
the  track  location  of  the  observations,  with  the  topographic 
features  in  the  background.  In  the  SAR  image,  features 
corresponding  to  bathymetric  variations  and  to  large-ampli¬ 
tude  internal  waves  can  be  seen  (Figures  3a  and  3b).  The 
bright  curve  marked  “bathymetry’'  in  Figure  3  corresponds 
to  topographic  contours  along  the  shoreline.  Hsu  el  ciL 
[2001]  has  noted  the  possibility  of  “seeing”  shallow  topo¬ 
graphic  features  using  SAR  images  owing  to  horizontal  tidal 
flow  convergence  over  rapidly  changing  topography.  The 
features  marked  B1  and  B2  are  curved  bright  lines,  which, 
in  light  of  the  model  simulations  to  be  presented,  possibly 
correspond  to  internal  hydraulic  jumps.  There  is  a  possibil- 


Depth  profile  from  (34N,  121 E)  to  (33.5Nr  123E) 


Figure  2.  Solid  line  extracted  topography  from  the  digital  atlas  of  Choi  [1999]  along  the  40-m  track. 
Dash-dotted  line  is  the  approximated  topography  with  an  analytical  function. 
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(a.) 


St 


S2 


S3 


(b.) 

Figure  3.  (a)  Location  of  features  in  RADARS  ATI  observations  on  8  August  1998.  B1  and  B2  are 

internal  bores.  SI  through  S4  indicate  a  propagating  solitary  wave  train.  A  bathymetric  feature  is  also 
present.  The  solid  black  line  represents  a  scale  of  46  km.  (b)  Zoom  of  internal  wave  solitary  wave  train  at 
SI,  S2,  and  S3.  The  solid  black  line  represents  a  scale  of  10  km. 


ity  that  a  better  pixel  resolution  might  show  that  the  jumps 
are  accompanied  by  a  few  waves.  The  bright  curved 
features  marked  SI  through  S4  are  solitary  wave  trains. 

[10]  The  solitary  wave  train  propagates  from  about  34°N, 
I22°E  to  33.4°N,  I24.5°E.  A  zoom  of  the  first  three 
resolved  solitary  wave  trains  is  shown  in  Figure  3b  (denoted 
SI,  S2,  and  S3).  The  spacing  between  the  bright  conver¬ 
gence  lines  increases  as  one  proceeds  from  SI  to  S3  (also 
through  S4,  although  this  is  difficult  to  see  from  Figure  3a). 
This  suggests  generation  in  the  coastal  shelfbreak  area. 
The  bright  and  dark  bands  mark  the  surface  convergence/ 
divergence  caused  by  the  current  fields  associated  with  the 
waves.  The  spacing  between  the  wave  packets  has  a  high 
degree  of  regularity  along  a  well-defined  direction.  The 
overall  concave  shape  of  the  wave  front  is  preserved 
along  the  travel  path,  which  again  suggests  that  these 
internal  wave  packets  originated  from  the  same  shelf 
break  location. 

[n  ]  As  the  solitary  wave  train  travels  away  from  the  point 
of  origin,  the  spacing  between  solitary  waves  in  the  train 
increases.  The  distance  between  the  first  and  the  sixth  wave 
packet  is  about  271  km.  Assuming  a  12.4-hour  semidiurnal 


tidal  period,  the  phase  speed  estimates  range  from  0.8  to 
l.l  m/s  along  the  track. 

2.3.  Tides 

[12]  The  dominant  barotropic  tidal  component  in  the 
Yellow  Sea  is  the  semidiurnal  M2  tide.  Pronounced  tidal 
variations  exist  along  the  Chinese  and  Korean  coasts,  where 
the  surface  elevation  changes  can  be  as  high  as  2  to  2.5  m 
(Figure  4)  [Naimie  et  ciL ,  2001],  There  are  strong  flood/ebb 
tidal  flow  points  around  the  Yangtze  River  on  the  Chinese 
coast  and  the  Han  River  on  the  Korean  coast  (Figure  la).  At 
these  points  the  tidal  flow  can  radiate  out  or  come  into  the 
region  [Yanagi  and  Takahashi ,  1993].  The  phase  of  the  flow 
at  these  locations  is  opposite.  These  two  regions  correspond 
to  areas  where  the  surface  elevation  contours  around 
amphidrous  points  meet  the  shore  [Yanagi  and  Takahashi , 
1993]. 

[13]  The  Shandong  area  of  our  studies  is  located  north  of 
the  Yangtze  River,  about  halfway  up  to  the  Shandong 
Peninsula.  There  the  shore  areas  still  have  pronounced  flood 
and  ebb  tides  [Yanagi  and  Takahashi ,  1993],  The  flood  tide 
component  shows  a  structure  that  is  indicative  of  a  dividing 
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Figure  4.  M2  surface  elevation  and  phase  structure,  in  cm/ 
deg,  simulated  with  a  barotropic  semidiurnal  tidal  model. 
The  phase  lines  originate  from  the  amphidrous  points  and 
are  color  coded  according  to  elevation  and  angle.  They  are 
color  coded  as:  25/0  (blue)  to  275/330  (red)  by  25/30  (cm/ 
deg).  Elevation  varies  from  25  cm  (blue)  to  275  cm  (red). 
Elevation  increment  is  25  cm  and  angle  is  30°.  The  lines  not 
connected  to  amphidrous  points  track  only  the  elevation. 


line  between  tidal  ellipses  that  are  cyclonic  in  the  northern 
part  of  the  Shandong  region  and  anti-cyclonic  in  the 
southern  pari  [Fang,  1994].  Along  the  shore  the  tidal 
magnitude  variability  of  the  area  ranges  from  1 .2  m/s  to 
0.3  m/s  [ Valencia ,  1987;  Fang ,  1994]. 

[u]  The  tidal  motion  across  the  topographic  variations 
generates  an  internal  tide  and  possibly  higher  frequency 
internal  waves,  which  propagate  away  from  the  generation 
site  [ Baines ,  1973,  1982].  These  waves  may  steepen 
through  nonlinear  effects,  after  which  frequency  and  ampli¬ 
tude  dispersion  sets  in  and  the  leading  edge  disintegrates 
into  high-frequency  solitary  wave  trains  [Lamb,  1994; 
Brandt  et  al ,  1997]. 

2.4.  Summer  Hydrography 

[i>]  For  the  Shandong  area,  hydrographic  surveys  are  not 
available.  As  a  result,  the  monthly  averaged  temperature 
and  salinity  distribution  from  the  National  Oceanographic 
Data  Center  (NODC)  database  was  used. 

[if*]  The  west  to  east  hydrographic  structure  at  35°N 
shows  variability  from  the  Chinese  to  the  Korean  coasts. 
Figure  5a  displays  some  temperature  change  along  the 
latitude  section  (121  °E  to  I24°E).  From  the  shallow  areas 
along  the  Chinese  coast  to  the  middle  of  the  Yellow  Sea  the 
thermocline  deepens,  from  about  10  m  to  25  m.  Climatolog¬ 
ical  6-month  averaged  salinity  distributions  contain  a  )0-m 
mixed  layer  and  exhibit  an  increase  in  salinity  from  longitude 
I21°E  to  124°E. 

[17]  At  the  outskirts  of  the  Shandong  study  area,  in  the 
vicinity  of  latitude  35°N  and  past  longitude  I24°E,  summer 
surveys  for  the  month  of  August  are  available.  Lie  [1986] 
has  extracted  temperature,  salinity,  and  sigma-t  distributions 


from  CTD  casts  made  in  August  of  1983  and  1984,  in  the 
southeastern  Hwanghae  (Yellow  Sea).  These  casts  give  us 
bounds  on  salinity,  temperature,  sigma-t,  and  the  general 
structure  at  the  outskirts  of  the  area.  From  the  NODC 
climatology  and  the  surveys  of  Lie  [1986],  generic  density 
structure  functions  are  constructed  that  have  the  form  of 

ct,  —  mz  +M[I  -  tanh(r  +  &)].  (I  ) 

where  z  is  depth  and  m,  A ,  and  b  are  arbitrary  constants 
selected  to  match  the  density  structure  of  the  region  for 
August. 

[is]  Figure  5b  shows  the  background  density  fields  that 
are  used  in  the  simulations.  Stratifications  with  pycnoclines 
at  depths  of  10,  15,  and  20  m  are  shown  (upper  40  m  only). 
For  these  stratifications,  solitary  waves  will  be  waves  of 
elevation/depression  in  the  shallow/deep  water,  depending 
on  the  pycnocline  location  [Liu  et  al .,  1998]. 

3.  Generation  and  Propagation  of  Solitary  Wave 

Trains 

3.1.  Model 

[19]  The  Lamb  [1994]  model  is  used  for  simulating  the 
generation  and  propagation  of  solitary  waves  in  the  Yellow 
Sea.  It  consists  of  the  Boussinesq  equations  with  the 
Coriolis  force  in  a  two-dimensional  cross-bank  plane.  The 
coordinates  are  defined  with:  x  as  across  bank,  y  as  along 
bank,  and  z  as  the  depth.  N ,  the  Brunt- Vaisala  frequency,  is  a 
function  of  depth  and  x  (range).  In  the  along-bank  direction, 
the  velocity  is  included  but  the  derivatives  are  neglected 
(2.5-dimensional  representation).  The  equations  of  the 
model  are 

V,  +  (V  •  V)V  -/V  x  k  =  -VP  -  k pg.  (2) 

p,  +  v  •  vp  =  0.  (3) 

v*v  =  0.  (4) 

where  \(u ,  v,  w)  is  the  velocity  vector,  V  is  the  three- 
dimensional  vector  gradient  operator,  subscript  t  denotes 
the  time  derivative,  p  is  the  density,  P  is  the  pressure,  g  is 
the  gravitational  acceleration,  /  is  the  Coriolis  parameter 
taken  as  9  x  10  5  s  1  for  a  latitude  of  40  degrees,  and  k  is 
the  unit  vector  along  the  z  direction.  In  the  three- 
dimensional  equations  (2)— (4)  the  partial  derivatives  with 

respect  toy  are  neglected,  i.e.,  —  (  )  -  0;  thus  equations  (2)- 

(4)  are  equivalent  to  equations  (la)-(ld)  in  Lamb's  [1994] 
paper. 

[20]  Before  the  equations  are  solved,  they  are  transformed 
to  a  terrain  following  coordinate  system  (sigma-coordinates) 
which  leads  to  higher  vertical  resolution  over  the  bank 
region.  The  equations  are  solved  over  a  domain  bounded 
below  by  the  topography  and  a  rigid  lid  above.  The  flow  is 
forced  by  specifying  a  semidiurnal  tidal  inflow  at  the  left 
boundary  of  the  form  VT  sin (^o/)  where  uj  is  the  M2  tidal 
frequency  assumed  to  have  a  1 2.4-hour  period.  The  strength 
of  the  tidal  current  in  the  shallow  water,  VT%  varies  between 
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Aug  Gen  Data  and  Clim  Temp  in  Yellow  Sea  at  35N  i-clim*124E  g-clim*122E.b-clirrv*121 
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Figure  5.  (a)  Temperature  profiles  from  climatology  at  35°N  and  I21°E  blue  line,  122°E  green  line,  and 
124°E  red  line,  (b)  Sigma-t  densities  with  pycnoclines  at  depths  of  10  m  (blue  line),  15  m  (red  line),  and 
20  m  (green  line). 


Yelbw  Sea  Shandong  August  Densities  A)st-g  osl-i  lOst-b 


0.3  and  1.2  m/s,  typical  of  values  in  the  Shandong  region 
(see  section  2.3).  On  the  right  boundary  an  outflow  condi¬ 
tion  is  used  [Lamb.  1994].  The  topographies  used  have 
already  been  discussed  in  section  2  (see  Figure  2).  The 
initial  conditions  are  a  state  of  rest  with  horizontally 
uniform  density  (section  2.4).  The  horizontal  grid  size  is 
20  m,  for  most  cases,  and  the  vertical  sigma-coordinate 


resolution  consists  of  60  grid  points.  The  time  step  is  tied  to 
the  Courant-Levi  condition  and  varies  from  I  s  to  5  s. 

3.2.  Dynamics  of  Generation  and  Propagation 

[21]  The  parameters  for  the  different  model  runs  are  given 
in  Table  I.  We  begin  by  considering  Case  2  in  detail.  For 
Case  2,  the  thermocline  depth  is  15  m,  the  peak  barotropic 
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Table  1.  Simulation  Parameters11 


Case 

h,  A  m 

//,  m 

F7,  m/s 

Topo 

dih  km 

1 

15 

70 

1.2 

a 

150 

T 

15 

70 

0.7 

a 

240 

3 

15 

70 

0.7 

a 

300 

4 

15 

70 

0.35 

a 

240 

5 

10 

70 

0.7 

a 

240 

6 

10 

40 

0.7 

a 

240 

7 

10 

40 

0.7 

b 

200 

X 

20 

70 

0.7 

a 

240 

9 

20 

70 

1.2 

a 

240 

'Hero  It,/  is  pycnocline  depth,  H  is  water  depth,  V,  is  the  tidal  strength, 
Topo  is  the  topography  type  (“a"  being  for  cases  with  a  finger,  and  “b” 
being  for  eases  with  the  finger  removed),  and  dtl  is  the  length  of  the 
computational  domain. 

tidal  velocity  is  0.7  m/s,  and  the  water  depth  is  70  in. 
Contour  plots  of  the  density  field  at  various  times  are  shown 
in  Figures  6a- 6e. 

[22]  Initially,  the  density  field  is  horizontally  homoge¬ 
neous  and  there  is  no  flow  in  the  .v  direction.  The  prescribed 
semidiurnal  tide,  on  the  left  boundary,  accelerates  the  fluid 
in  the  off-bank  direction.  The  presence  of  a  ridge  in  the 
topography  results  in  the  generation  of  two  initial  depres¬ 
sions  during  the  off-shelf  flow.  Each  of  these  splits  into  two 
waves,  each  propagating  in  opposite  directions.  In  Figure  6a 
the  small  depression  over  the  ridge  is  a  combination  of  two 
of  these,  and  is  predominantly  comprised  of  the  off-shelf 
depression  formed  over  the  first  drop  in  the  topography.  The 
shoreward  propagating  depression  formed  over  the  first 
drop  has  steepened  to  form  a  hydraulic  jump.  Its  shoreward 
propagation  has  been  halted  for  a  time  by  the  supercritical 
tidal  flow  on  the  shelf.  During  on-shelf  tidal  flow  the 
bottom  of  the  depression  is  pushed  up  and  a  large  jet  of 
water  rises  well  above  the  bottom.  In  the  end,  a  wave  of 
elevation  is  formed.  One  such  wave  is  formed  each  tidal 
period  as  can  be  seen  in  Figure  6b,  which  shows  the  density 
field  at  42  hours  (about  3.4  tidal  periods).  Here  the  leading 
small  depression  at  about  125  km  is  the  off-shelf  propagat¬ 
ing  depression  formed  to  the  lee  of  the  ridge  during  the  first 
tidal  period.  The  first  wave  train  to  come  off  the  shelf,  at 
about  110  km,  was  formed  from  the  off-shelf  propagating 
depression  formed  over  the  first  drop,  between  the  shelf  and 
the  ridge,  during  the  first  tidal  period.  During  the  second 
period  of  off-shelf  flow,  it  was  just  to  the  right  of  the  ridge  at 
the  location  where  the  second  initial  depression  was  formed 
in  the  first  tidal  period.  Hence  the  creation  of  a  depression  to 
the  lee  of  the  ridge  can  only  be  seen  during  the  first  tidal 
period. 

[23]  After  four  tidal  cycles  a  fourth  solitary  wave  train  is 
about  to  form  from  a  large  hydraulic  jump  at  around  30  km 
(Figure  6c).  The  second  solitary  wave  train,  from  the  left,  in 
Figure  6b  has  moved  to  100  km  (Figure  6c).  Note  that  the 
spacing  between  the  waves  in  the  first  packet  has  increased. 
This  is  the  first  train  to  come  off  the  shelf  and  has  evolved  to 
a  structure  in  which  the  individual  waves  in  each  packet  get 
farther  a  pan. 

[24]  Each  tidal  cycle  generates  wave  propagating  on  the 
shelf  and  a  wave  packet  propagating  away  from  the  shelf. 
This  behavior  is  seen  in  other  areas  [ApeL  2003;  Brandt  et 
a /.,  1997;  Warn-Varnas  et  ui,  2003].  The  waves  of  elevation 
on  the  shelf  travel  toward  the  shore  and  oscillate  with  the 
semidiurnal  tidal  motion.  The  amplitudes  of  the  waves 


decrease  through  dynamical  evolution,  and  they  eventually 
disappear.  At  63  hours  or  5.1  semidiurnal  tidal  cycles  into 
the  simulation,  there  are  four  well-developed  wave  packets, 
with  a  fifth  starting  to  form  (Figure  6d).  Note  that  the  first 
three  wave  packets  from  the  shelf  show  a  dramatic  increase 
in  amplitude.  This  is  largely  due  to  the  response  over  the 
shelf  edge  increasing  in  time,  as  discussed  by  Lamb  [1994]. 
The  third  and  fourth  packets  from  the  shelf  are  more  similar 
in  size  (when  compared  at  similar  stages  of  their  evolution). 
The  individual  waves  in  each  packet  grow  in  size  for  a 
while  and  then  start  to  decay.  They  also  get  farther  apart. 
This  is  particularly  apparent  in  the  first  two  packets  away 
from  the  shelf  shown  in  Figure  6e.  The  decrease  in 
amplitude  is  a  physical  phenomenon.  Numerical  dissipation 
and  boundary  effects  are  minor  influences. 

[25]  Note  the  recovery  of  the  isopycnals  toward  the 
original  equilibrium  position  before  the  next  depressions 
of  the  new  solitary  wave  train  in  the  following  tidal  cycle. 
At  the  end  of  one  tidal  period,  each  isopycnal  must  have  the 
same  average  height  as  it  had  at  the  start  of  the  simulation. 
In  the  case  of  waves  of  depression  propagating  away  from 
the  generation  site,  conservation  of  mass  implies  that  the 
isopycnals  must  lie  above  their  initial  height  somewhere  in 
the  domain  after  a  full  tidal  period  (this  is  not  true,  for 
example,  half  a  tidal  period  after  the  start  of  off-shelf  flow). 
This  occurs  at  the  generation  site  where  the  dynamics  cause 
the  fluid  to  rebound  and  overshoot  its  initial  height  after 
being  forced  downward  during  off-shelf  flow  [Lamb,  1994]. 
The  spacing  between  the  wave  trains  tends  to  decrease 
toward  the  right  boundary,  indicating  that  the  phase  velocity 
decreases  as  the  trains  propagate  to  the  far  right.  This  is 
partly  a  consequence  of  the  fact  that  the  wave  propagation 
speed  is  an  increasing  function  of  amplitude. 

[26]  The  model  simulations  establish  the  generation  of 
solitary  waves  in  the  shelfbreak  area  and  the  off-shelf 
propagation  of  solitary  wave  packets.  The  solitary  waves 
are  generated  by  semidiurnal  tidal  motion  over  topographic 
changes  at  the  shelfbreak.  SAR  imagery  also  indicates 
solitary  wave  packets  propagating  off  the  shelf  (sections 

2.2  and  3.6),  suggesting  the  generation  of  the  wave  packets 
at  the  shelfbreak.  Solitary  wave  packet  generation  and 
propagation  can  occur  in  other  areas  of  the  Yellow  Sea. 
Hsu  et  al.  [2000]  tracked  solitary  wave  packets  propagating 
from  east  to  west.  Most  of  these  solitary  wave  packets 
indicate  an  origin  southwest  off  the  tip  of  Korea.  Some  of 
the  solitary  wave  packets  have  an  origin  in  the  coastal  area 
off  China. 

[27]  The  strength  of  the  semidiurnal  tide,  together  with 
the  structure  of  topography,  influences  the  size  of  the  initial 
depressions  formed  over  the  bank  edge  during  the  off-bank 
flow.  Lamb  [1994]  estimated  the  size  of  the  initial  pycno¬ 
cline  depressions  that  arise  from  variation  in  the  prescribed 
tidal  strength,  V-f,  and  the  topographic  structure.  He  let  the 
depression  of  volume  V  formed  over  the  shelf  slope  have  a 
characteristic  length  L  and  depth  D.  The  length  of  the 
depression  is  determined  by  the  width  of  the  shelf  or  bank 
edge  and  an  analysis  of  the  volume  displacement  caused  by 
the  tide  showed  that  the  amplitude  of  the  initial  depression 
scales  as  D  Vt!L. 

[2s]  The  tidal  parameter,  VT,  was  varied  from  0.35  m/s  to 

1.2  m/s.  These  values  reflect  possible  conditions  in  the 
Shandong  area  (see  section  2.3).  The  consequences  of 
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Figure  6.  Model  predicted  evolution  of  solitary  wave  trains  over  semidiurnal  tidal  periods.  Sigma-t 
density  field  for  Case  2  in  Table  I  at  semidiurnal  tidal  periods  of  (a)  0.48,  (b)  3.4,  (c)  4.0,  (d)  5.1,  and 
(e)  6. I . 


varying  the  tidal  strength  can  be  illustrated  by  comparing 
Cases  2  and  4  in  Table  1.  The  only  difference  between  the 
two  cases  is  the  strength  of  the  tidal  forcing  which  decreases 
from  VT  =  0.7  m/s  for  Case  2  to  VT  =0.35  m/s  in  Case  4.  A 
comparison  of  the  results  at  42  hours  is  shown  in  Figure  7. 
For  a  tidal  forcing  of  VT-  0.7  m/s  the  first  two  wave  trains 
to  come  off  the  shelf  are  located  near  105  km  and  78  km. 
The  second  wave  train  has  a  leading  amplitude  of  about  7  m 
and  a  wavelength  of  475  m  (distance  between  wave  crests) 
near  the  front  of  the  train  (Figure  7a).  For  a  tidal  forcing  of 
VT  =0.35  m/s  the  previous  first  wave  train  to  come  off  the 
shelf  does  not  exist.  The  second  solitary  wave  train  has  a 


leading  amplitude  of  around  4  m  and  a  wavelength  of  400  m 
at  a  distance  of  around  78  km  (Figure  7b).  As  the  tidal 
forcing  increases,  the  volume  of  water  displaced  increases 
and  results  in  larger  initial  depressions.  As  a  result,  the 
amplitudes  and  wavelengths  of  the  initial  solitary  wave  train 
that  develops  are  larger.  For  the  second  solitary  wave  train, 
there  is  an  increase  in  amplitude  as  the  tidal  forcing 
increases.  The  wavelength  exhibits  a  small  increase  (from 
400  m  to  475  m)  with  increasing  tidal  forcing  from  0.35  m/s 
to  0.7  m/s. 

[2<>]  For  a  thermocline  located  at  20  m.  Cases  8  and  9  in 
Table  1,  a  corresponding  analysis  was  conducted  for  tidal 
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Figure  7.  Simulations  at  3.22  semidiurnal  tidal 
periods  with  tidal  forcing  of  (a)  Vr  ~  0.7  m/s  and  (b)  VT  = 
0.35  m/s. 


forcings  of  0.7  m/s  and  1.2  m/s.  The  results,  for  the  second 
wave  train  to  come  off  the  shelf  at  around  83  km,  yield 
amplitudes  of  8  m  and  wavelengths  of  480  m  for  both  tidal 
forcings.  These  are  similar  values  to  the  previous  case  for  a 
tidal  forcing  of  0.7  m/s.  This  suggests  that  once  the  tidal 
strength  reaches  a  certain  magnitude  or  saturation  value  (as 
for  the  cases  of  Vr  =  0.7  m/s  and  Vr  =  1.2  m/s),  the 
amplitudes  and  wavelengths  of  the  second  solitary  wave 
train  (and  higher)  reaches  a  saturation  limit  also,  and  they 
tend  not  to  increase  further. 

|>]  Near  the  saturation  value  of  amplitude  and  wave¬ 
length,  the  most  pronounced  effect  of  tidal  forcing  is  on 
the  first  wave  train  that  comes  out.  For  a  forcing  of  Vr  = 
0.35  m/s,  there  is  no  first  solitary  wave  train.  For  forcings 
of  Vr  -  0.7  m/s  and  VT  =  1.2  m/s  the  first  wave  train 


increases  in  amplitude  in  relation  to  the  tidal  forcing 
magnitude. 

3.3.  Effects  of  Topographic  Variations 

[3 1 ]  The  topographic  structure  in  the  Shandong  area  is 
described  in  section  2.1 .  Two  transects  across  the  shelf  slope 
have  been  considered.  Both  consist  of  an  initial  drop,  a  rise, 
and  then  a  second  drop  to  their  final  depths  of  40  and  70  m. 
Results  from  a  simulation  using  a  final  depth  of  40  m  (Case 
6  in  Table  I )  are  shown  in  Figure  8a  at  42  hours.  This  can  be 
compared  with  Figure  6b,  which  shows  results  of  Case  2  at 
the  same  time.  Case  6  uses  a  thermocline  depth  of  10  m 
because  a  thermocline  depth  of  15  m  did  not  result  in  the 
formation  of  a  solitary  wave  train.  The  position  of  the  first 
wave  train  to  come  off  the  shelf,  for  the  70  m  depth,  in 
Figure  6b  is  just  past  100  km.  For  the  40-m-depth  case  the 
first  wave  train  to  come  off  the  shelf  is  at  about  95  km, 
suggesting  a  smaller  phase  speed.  The  next  wave  train  has  a 
double  structure,  with  two  solitary  wave  train  staictures 
close  together  at  around  70  km  (Figure  8a).  The  first  jump 
(and  solitary  wave  train  to  appear)  evolved  from  the  right 
propagating  depression  in  the  lee  of  the  first  drop  reinforced 
by  the  left  propagating  depression  from  the  second  rise. 
Subsequently,  for  the  following  tidal  cycles  the  modulation 
of  the  right  propagating  depression  (from  the  first  drop)  by 
the  left  propagating  depression  from  the  second  rise  was 
such  that  a  double  structure  solitary  wave  train  resulted.  In 
Case  2,  with  a  deep  water  depth  of  70  m,  the  two  off-shelf 
propagating  depressions  were  superimposed  on  each  other, 
resulting  in  a  single  jump  and  wave  train  formed  during 
each  tidal  period.  A  second  double  structured  solitary  wave 
train  is  starting  to  fonn  at  about  35  km.  The  wave  trains  that 
have  evolved  in  the  40-m  channel  contain  fewer  waves  and 
are  much  smaller  than  the  trains  in  the  70-m  channel  shown 
in  Figure  6b.  On  the  other  hand,  on  the  shelf  the  internal 
bore  structures  are  more  pronounced  in  the  40-m  channel 
simulations.  This  is  due  to  the  pycnocline  being  centered 
above  the  top  of  the  shelf,  resulting  in  more  stratification  in 
the  shallow  water  in  Case  6  than  in  Case  2. 

[32]  For  the  40-m  channel  the  effects  of  the  second 
topographic  rise  (or  finger)  were  investigated  by  removing 
it.  All  other  aspects  of  the  simulation  were  unchanged 
(Case  7). 

[33]  The  results  at  hour  42  are  shown  in  Figure  8b.  The 
double  structured  wave  trains  have  disappeared  as  expected. 
In  Case  6  the  first  wave  train  is  at  95  km  at  t  =  42  hours. 
This  wave  train  is  not  present  in  Case  7,  indicative  of  the 
fact  that  the  left  propagating  depression  from  the  second  rise 
exerted  a  reinforcement  effect. 

[34]  Along  the  SAR  imagery  and  simulation  track,  there 
is  a  depth  gradient  from  shallow  to  deep  water,  not  only  in 
the  internal  solitary  wave  propagation  direction,  but  also  in 
the  along  crest  direction  with  the  southern  end  of  the  train 
being  over  shallower  waters  than  the  northern  end.  As  the 
internal  solitary  waves  enter  the  shallow  waters  (depths  of 
around  40  m)  from  the  generation  point,  they  disappear.  The 
disappearance  of  internal  solitary  waves  has  been  observed 
in  regions  where  the  waves  shoal  into  a  coastal  area.  In  such 
situations,  the  pycnocline  is  at  about  mid-depth  [Porter  and 
Thompson ,  1999;  Porter  et  al. ,  2001].  This  is  not  the  case 
here,  since  the  disappearance  of  internal  waves  is  near  the 
generation  area.  Our  model  simulations,  however,  indicate 
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(b) 

Figure  8.  Simulations  at  3.22  semidiurnal  tidal  periods 
with  two  topographies:  (a)  Case  6  in  Table  1  and  (b)  Case  7 
in  Table  1 . 


that  for  a  pycnocline  depth  of  20  m  in  a  40-m  channel, 
internal  solitary  waves  do  not  evolve  as  predicted  by  weakly 
nonlinear  theory  for  symmetric  stratifications  [Liu  el  aL, 
1998].  For  other  pycnocline  depths  in  40-m  depth  channels, 
such  as  10  m,  solitary  waves  do  evolve.  Such  conditions  are 
possible  during  some  portion  of  the  summer  months. 

3.4.  Dispersion  Variations  in  the  70-m-Depth  Channel 

[35]  The  approach  for  the  analysis  of  the  wavelength  and 
amplitude  of  individual  wave  trains  is  illustrated  by  the 
schematic  diagram  shown  in  Figure  9.  Waves  of  depression 
have  peak  horizontal  surface  currents  in  the  center  of  the 
wave  in  the  direction  of  propagation.  Thus  peak  surface 
convergence  (max  duidx)  occurs  in  the  front  half  of  the 
wave  while  peak  surface  divergence  occurs  in  the  second 
half  of  the  wave.  The  distance  between  the  extrema  in  Du! Ox 


is  called  D.  Regions  of  strong  convergence  and  divergence 
can  appear  as  bright  and  dark  lines  in  SAR  observations. 
The  wavelength,  X,  is  defined  as  the  distance  between 
successive  centers  of  troughs.  The  wavelength  will  be  close 
to  the  distance  between  the  convergences  or  bright  lines  in 
SAR  observations.  For  well-separated  solitary  waves  it  is 
the  distance  between  the  waves,  not  a  measure  of  the  length 
of  an  individual  wave.  Except  for  very  large  times,  the 
waves  are  not  well  separated  in  the  model  results.  This 
definition  of  the  wavelength  is  used  in  order  to  facilitate  the 
interpretation  of  the  Fourier  spectra  and  the  relation  of  the 
model  simulations  to  SAR  observations.  Wave  amplitudes 
are  defined  as  the  vertical  displacement  of  an  isopycnal 
from  just  ahead  of  an  individual  wave.  An  isopycnal  in  the 
middle  of  the  thermocline  is  used,  in  this  case  with  a  value 
of  sigma-t  =  22.5.  This  definition  is  somewhat  problematic, 
particularly  for  the  leading  wave,  as  the  wave  train  is  in  the 
depression  of  the  internal  tide  and  the  internal  tide  changes 
shape  with  time.  In  particular,  the  isopycnal  depressions 
rise,  after  the  solitary  wave  packet,  and  then  return  to 
equilibrium  (e.g.,  see  Figure  6).  The  dispersive  behavior 
of  the  second  solitary  wave  train  (to  come  off  the  shelf)  for 
Case  2  in  Table  1  has  been  analyzed  at  30,  42,  55,  and 
67  hours.  This  is  the  case  with  a  tidal  forcing  of  VT=  0.7  m/s 
and  a  pycnocline  depth  of  hj  =  15  m  depicted  in  Figure  6. 
Figure  10a  shows  the  evolution  of  the  solitary  wave  trains 
from  30  to  55  hours.  At  42  hours  the  second  train  is  located 
near  80  km.  At  55  hours  it  has  moved  to  near  120  km  and  has 
noticeably  changed  its  structure.  An  enlargement  of  the 
second  solitary  wave  train  is  displayed  in  Figure  10b.  A 
distinct  change  in  the  structure  of  this  train  can  be  seen  for  the 
time  interval  of  42  to  55  hours.  At  55  hours  the  back  of  the 


Figure  9.  Relationships  of  solitary  waves  to  surface 
signatures.  Arrow  shows  propagation  direction.  Upper 
curve  is  a  typical  isopycnal.  Lower  curve  is  duidx  where 
u  is  the  surface  velocity  (waves  are  propagating  to  the  right). 
Positive/negative  values  are  regions  of  surface  convergence/ 
divergence.  D  is  the  distance  between  the  maxima  in  surface 
convergence/divergence  which  appear  as  bright  and  dark 
lines  in  SAR  observations.  The  distance  between  waves,  X, 
will  be  referred  to  as  the  wavelength.  In  SAR  images  it  is 
the  distance  between  bright  lines. 
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Figure  10.  (a)  Sigma-t  =  22.5  isopycnal  at  times  of  30,  42,  55,  and  67  hours,  for  Case  2  in  Table  1. 
(b)  Sigma-t  -  22.5  expanded  view  of  isopycnals  at  times  of  30,  42,  55,  and  67  hours,  (c)  Amplitude- 
wavelength  dispersion  diagram  for  second  solitary  wave  train,  showing  wavelength  versus  amplitude, 
for  Case  2  in  Table  1  at,  30,  42,  55,  and  67  hours.  Solid  lines  are  polynomial  fits. 
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Dispersion  for  3  cases  of  thermocline  depth 


Figure  11.  Dispersion  diagrams  at  50  hours  exhibiting  effects  of  pycnocline  variations  on  second 
solitary  wave  train,  showing  wavelength  versus  amplitude,  pycnocline  at  10-m  curve  c,  1 5-m  curve  b, 
and  20-m  curve  a. 


train  has  larger  wavelengths  than  the  front  of  the  train.  This 
trend  is  even  more  pronounced  and  dominant  at  67  hours.  The 
wavelengths  as  a  function  of  amplitude  distributions  of  the 
solitary  wave  train  at  the  four  times  are  shown  in  Figure  10c. 
At  30  hours  the  amplitudes  vary  from  about  9.5  to  0.2  m  and 
the  wavelengths  decrease  from  about  250  to  200  m.  At 
42  hours  the  amplitudes  range  is  slightly  less  and  the  wave¬ 
lengths  have  increased.  They  now  vary'  from  425  at  the  front 
of  the  wave  train  to  350  m  at  the  back.  The  large  decrease  in 
the  amplitude  of  the  leading  wave  is  due  to  a  drop  in  the 
size  of  the  jump  at  the  front  of  the  wave.  At  55  hours  the 
structure  of  the  second  wave  train  changes,  in  the  sense  that 
the  wavelengths  tend  to  increase  from  the  front  to  the  back  of 
the  wave  train  with  the  amplitudes  decreasing  as  before.  To 
distinguish  between  different  behaviors,  a  wave  train  with 
wavelengths  decreasing  from  front  to  back  will  be  called  a 
type  A  configuration.  Similarly,  wave  trains  whose  wave¬ 
lengths  increase  from  front  to  back  will  be  called  a  type  B 
configuration.  The  wavelength  variation  is  from  550  to 
1000  m.  At  67  hours  the  amplitudes  decrease,  as  before,  and 
the  wavelengths  show  more  increase  than  at  55  hours  from 
front  to  back,  ranging  from  625  to  2300  m.  At  the  back  of 
the  wave  train  the  distances  between  the  centers  of  the 
depression  increase  appreciably  as  evidenced  in  Figure  6d. 

[36]  Away  from  the  shelf  and  topographic  changes,  con¬ 
figuration  B  appears.  When  the  configuration  B  occurs, 
the  phase  speeds  of  the  depressions  in  the  rear  of  the 
train  decrease.  Then  the  spaces  between  the  depressions, 
in  the  rear,  increase  and  the  amplitude  of  the  depressions 


decreases.  As  the  amplitude  decreases,  the  wavelength 
increases. 

[37]  In  the  second  wave  train  (to  come  off  the  shelf),  the 
location  of  the  center  of  the  first  depression  was  tracked  in 
time.  A  straight  line  fit  yields  a  propagation  speed  of 0.83  m/s. 
This  is  consistent  with  estimates  from  SAR  observations 
(section  3.6),  which  yield  values  of  0.8  to  1.1  m/s.  The 
model  results  were  obtained  with  the  summer  stratification 
derived  from  climatology  and  available  data  (section  2.4). 
A  thermocline  depth  of  15  m.  Case  2  in  Table  1,  was 
found  to  yield  the  closest  results  to  the  SAR  measurements 
(sections  3.6  and  3.7). 

3.5.  Effects  of  Pycnocline  Variation 

[3«s]  The  summer  hydrography  (section  2.4)  was  used  to 
determine  the  density  field  used  in  the  model  simulations. 
The  relevant  parameter  that  describes  the  density  structure 
for  typical  summer  conditions  is  the  pycnocline  depth,  ht). 

[39]  For  a  tidal  forcing  of  Vt~  0.7  m/s,  we  considered  the 
sensitivity  of  the  model  results  to  the  stratification  by 
comparing  the  results  for  the  following  thermocline  depths: 
10  m,  15  m,  and  20  m  (Cases  5,  2,  and  8,  respectively,  in 
Table  l).  The  wavelength  versus  amplitude  dispersion 
analysis  of  the  second  solitary  wave  train  (to  come  off  the 
shelf)  at  50  hours  is  shown  in  Figure  11,  for  the  three 
thermocline  locations.  For  a  20-m  location.  Case  8  in 
Table  I,  the  amplitudes  range  from  5  m  to  1.5  m  and  the 
wavelengths  from  320  m  to  400  m  (curve  in  Figure  1 1). 
The  wavelengths  increase  from  the  front  to  the  back  of  the 
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Table  2.  Parameters  of  SAR  Observations 


Bore  1 

Bore  2 

Packet  1 

Packet  2 

Packet  3 

Packet  4 

Packet  5 

Distance,  km 

0 

36  46 

72  93 

132 

178 

226 

272 

Phase  speed  C,  m/s 

0.8  1.03 

0.8  1.05 

1 .098 

1.03 

1 .075 

1.03 

Dominant  wavelength,  m 

630 

930 

1600 

2300 

solitary  wave  train,  indicating  the  presence  of  type  B 
configuration.  For  a  15-m  thermocline  location.  Case  2  in 
Table  1,  the  wavelengths  increase  from  530  m  to  700  m  and 
the  amplitudes  vary  from  9.5  m  to  1.5  m  (curve  “b”  in 
Figure  II).  The  wavelength  trend  also  exhibits  a  type  B 
behavior.  For  a  10-m  thermocline  location,  Case  5  in  Table 
l,  the  wavelength  data  exhibit  scatter  about  the  fitted  curve 
that  show  wavelengths  of  around  525  m  in  front  of  the  train, 
510  m  in  the  back,  and  570  m  in  the  middle,  suggesting  a 
transition  between  types  A  and  B  behavior.  The  amplitudes 
vary  from  12  m  to  about  4.0  m  (curve  tkc”  in  Figure  1 1). 
Overall,  the  largest  wavelengths  occur  for  Case  2,  the  15-m 
thermocline,  shown  as  curve  “b”  in  Figure  11.  This 
suggests  a  15-m  optimum  thermocline  location  for  maxi¬ 
mum  simulated  solitary  wave  intensity.  The  phase  speeds 
vary  from  0.73  m/s  to  0.82  m/s.  The  15-m  thermocline 
location  yields  a  phase  speed  of  0.82  m/s,  which  is  a  lower 
bound  in  relation  to  SAR  observations  (Table  2). 

3.6.  Comparisons  With  SAR  Observations  and 
Discussion 

[40]  The  SAR  observations  indicate  a  wave  packet 
propagating  away  from  shore  along  the  track  shown  in 
Figure  la.  In  the  direction  of  the  travel  path,  there  is  a 
depth  gradient,  from  shallow  to  deep  water,  not  only  in  the 
propagation  direction  but  also  in  the  along-crest  direction 
with  the  southern  end  of  the  solitary  wave  train  being  over 
shallower  water  than  the  northern  end.  In  addition,  as  the 
south  end  of  the  wave  crest  enters  shallow  waters,  its 
signature  tends  to  disappear  (M.  Teixeira,  personal  com¬ 
munication,  2002).  The  southern  end  of  the  signature  is 
located  around  the  40-m  depth  region.  Model  simulations 
indicate  that  for  a  15-m  pycnocline  in  a  water  depth  of 
40  m,  no  solitary  waves  develop  (section  3.3).  As  a  result, 
we  will  only  compare  the  70-m  channel  simulations  with 
the  SAR  data. 

[41]  Along  the  track  there  are  the  signatures  of  two 
internal  hydraulic  jumps  and  five  solitary  wave  trains, 
shown  in  Figure  3a.  A  zoom  of  the  first  resolved  solitary 
wave  train  is  shown  in  Figure  3b,  on  a  20  km  by  20  km  grid. 
The  distances  between  the  bright  lines  represent  the  spacing 
between  the  train  convergence  locations  (see  Figure  9),  and 
indicate  the  wavelength  between  the  centers  of  the  troughs 
associated  with  the  depression  of  the  isopycnals. 

[42]  A  relationship  between  the  half-width  L  of  a  solitary 
wave  and  the  distance  D  between  maximum  points  of 
surface  convergence  and  divergence  in  a  solitary  wave 
can  be  derived.  The  derivation  is  based  on  the  hyperbolic 
secant  representation  of  a  solitary  wave  [Small  et  a /.,  1999] 
and  gives  the  relationship 

L  =  D/1.32.  (5) 

We  have  estimated  the  distance  D  from  the  RADARSAT1 
ScanS A R  observations  for  the  first  wave  train  shown  in 


Figure  3b.  Assuming  a  resolution  of  100  m  per  pixel  and 
measuring  the  first  wavelength  in  the  train,  one  obtains  L  = 
75.7  m  (see  Texeira  et  al.  [2005]  for  more  details).  The  error 
is  large,  being  ±37  m  (assuming  an  error  of  50  m  in  the  SAR 
image). 

[42]  The  results  of  two-dimensional  Fourier  spectral  anal¬ 
ysis  of  the  SAR  images  are  summarized  in  Table  2,  where 
the  solitary  wave  trains  are  labeled  from  left  to  right.  The 
listed  wavelengths  are  for  the  most  intense  spectral  peaks 
that  occur  for  packets  1,  2,  3,  and  4.  For  packet  5,  there  is 
not  enough  signal  above  background.  The  Fourier  analysis 
reflects  the  wavelength  around  which  most  of  the  energy  is 
concentrated.  The  set  of  wavelengths  that  we  obtained  for 
the  solitary  wave  trains  is  630  m,  930  m,  1 600  m,  and  2300  m 
(Table  2). 

[44]  The  locations  of  the  internal  bores  and  solitary  wave 
packets  are  listed  in  Table  2  relative  to  the  first  bright  spot 
appearance  or  beginning.  Bore  1  is  located  on  the  shelf- 
break.  It  is  assumed  that  each  solitary  wave  packet  is 
generated  at  intervals  of  the  semi-diurnal  tidal  period  or 
12.4  hours.  With  this  assumption,  the  phase  speeds  are 
estimated  along  the  track.  They  range  from  0.8  to  1 .098  m/s. 

[45]  For  comparison  with  SAR,  the  simulation  with  a  1 5-m 
thermocline.  Case  2  in  Table  1 ,  is  used.  This  is  the  simulation 
whose  second  solitary  wave  train  (that  comes  off  the  shelf) 
has  the  largest  wavelengths  at  four  tidal  cycles  (see  Figure  1 1). 
The  wave  train  is  located  at  100  km  (Figure  6c).  The 
calculated  spectra  are  shown  in  Table  3.  The  dominant 
wavelength  in  the  vicinity  of  100  ktn  is  around  640  m  and 
is  practically  at  the  630-m  value  contained  in  the  data 
(Tables  2  and  3).  The  agreement,  however,  changes  as 
more  and  more  solitary  wave  trains  are  generated  and 
move  through  the  100  km  location.  Table  3  shows  the 
calculated  wavelengths  at  the  various  horizontal  locations 
(labeled  from  left  to  right)  for  different  trains  that  are 
generated  during  semidiurnal  tidal  cycles.  The  first  train 
that  propagates  away  from  the  shelf  break  region  behaves 
differently  from  the  other  trains  and  represents  a  different 
dynamical  state  of  evolution.  Consider  this  to  be  the 
initial  adjustment  state. 


Table  3.  Model  Results  for  Case  2  in  Tabic  I  (240  km)  for  f  irst 
Through  Fifth  Wave  Trains*1 


1 

2 

3 

4 

5 

Tidal  cycle 

4 

4 

4 

Distance,  km 

65 

100 

130 

Wavelength,  m 

330 

640 

1300 

Tidal  cycle 

5.1 

5.1 

5.1 

5.1 

Distance,  km 

60 

100 

140 

165 

Wavelength,  in 

270 

490 

640 

2700 

Tidal  cycle 

6.1 

6.1 

6.1 

6.1 

6.1 

Distance,  km 

63 

102 

143 

178 

195 

Wavelength,  m 

335 

420 

810 

640/2300 

3300 

:iThe  slash  indicates  another  wavelength  that  is  present.  Distance  is 
measured  from  the  sheltbreak. 
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Table  4.  Model  Results  for  Case  3  in  at  Tidal  Cycle  7.1  for  First  Through  Sixth  Wave  Trains* 


1 

2 

3 

4 

5 

6 

Distance,  km 

95 

138 

179 

211 

224 

Wavelength,  m 

340 

640 

1050/2100 

550/2300 

1300/2500 

‘The  slash  indicates  another  wavelength  that  is  present.  Distance  is  measured  from  the  shelfbreak. 


[46]  The  solitary  wave  trains  propagate  away  from  the 
shelf  through  the  various  horizontal  locations  or  check 
points  at  which  the  spectra  are  calculated.  Similar  loca¬ 
tions  are  chosen  for  comparing  model  predicted  wave¬ 
lengths  with  SAR  measurements.  At  around  100  km  the 
model  predicted  wavelength  varies  from  640  m  to  420  m 
(disregarding  the  adjustment  state)  as  the  tidal  cycle 
progresses  from  4  to  6.1  (Figure  6e).  The  measurements 
(Table  2)  show  a  630-m  wavelength  at  the  location.  In 
the  vicinity  of  130  km  to  140  km  the  model  results  yield 
a  variation  of  640  m  to  810  m  while  the  data  show  930  m.  At 
around  170  km  the  modeled  wave  train  is  displaying  an 
increased  spacing  between  waves  that  is  most  pronounced 
toward  the  back  of  the  wave  packet.  The  resultant  wavelength 
due  to  this  increased  spacing  is  around  2300  m  (Table  3).  The 
corresponding  wavelength  in  the  measurements  is  around 
1600  m  (Table  2).  This  is  a  situation  where  the  model  over¬ 
predicts  the  wavelength  instead  of  underpredicting  it.  This 
suggests  a  change  in  the  behavior  of  the  measured  solitary 
wave  trains  from  type  A  to  type  B  configuration  that  results  in 
a  sudden  increase  in  wavelength  size.  The  model  results 
at  ranges  greater  than  170  km  also  indicate  such  a 
possibility. 

[47]  The  change  in  behavior  of  solitary  wave  trains  is 
investigated  further  by  considering  a  larger  horizontal 
domain  in  which  the  solitary  wave  trains  can  propagate 
further  relative  to  the  previous  simulation.  Case  3,  in  Table  1 . 
The  horizontal  domain  is  increased  from  240  km  to  300  km. 
The  resolution  is  60  grid  points  in  the  vertical,  12,000  grid 
points  in  the  horizontal,  for  both  cases.  The  depth  of  the  ocean 
is  70  m  for  both  cases.  The  horizontal  grid  spacing  changes 
from  20  m  to  25  m  for  the  300-km  horizontal  distance,  Case  2 
in  Table  1 . 

[48]  Table  4  shows  the  analyzed  Fourier  wavelengths  for 
this  case  (Case  3  in  Table  1).  The  simulations  of  Table  4 
extend  to  7.1  tidal  cycles.  The  dominant  large  wavelengths 
show  an  abrupt  increase  from  wave  trains  1  and  2  to  wave 
trains  3  and  4.  This  abrupt  increase  is  also  present  in  the 
data  (Table  2).  This  reinforces  the  hypothesis  of  two 
behavior  states,  A  and  B. 

[49]  The  variation  of  modeled  wavelengths  is  due  to 
several  factors.  As  the  solitary  wave  trains  travel  through 
the  region,  the  isopycnals  tend  to  return  to  an  equilibrium 
position  that  becomes  progressively  higher  than  the  original 
(Figure  6). 

[50]  Another  source  of  discrepancy  is  the  three  dimen¬ 
sionality  of  the  problem.  The  curvature  effects  in  the 
2.5-dimensional  model  that  neglects  the  derivatives  per¬ 
pendicular  to  the  2-D  plane  are  not  represented,  as  well 
as  any  cross-plane  dynamical  processes  that  involve 
transport  and  exchange.  The  neglect  of  three-dimensional 
effects  probably  plays  a  role  in  the  different  character¬ 
istics  of  the  solitary  wave  trains  that  propagate  off  the 
shelf  break. 


[51]  Another  parameter  to  consider  is  horizontal  grid 
resolution.  The  results  shown  in  Tables  3  and  4  were 
obtained  with  horizontal  resolution  of  20  m  and  25  m, 
respectively.  The  20-  to  25-m  horizontal  resolutions  resolve 
the  internal  solitary  wave  wavelengths  that  range  from  630  m 
to  2300  m  in  the  SAR  observations  (Table  2).  The  20-  to  25-m 
grid  horizontal  intervals  are  also  within  the  100-m  pixel 
resolution  of  the  SAR  data  and  the  half-width,  L. 

[52]  The  model  simulations  exhibit  the  presence  of  two 
behavior  configurations,  A  and  B,  that  have  corresponding 
solitary  wave  train  characteristics  These  states  of  behavior 
are  evolved  by  the  dynamics  of  the  model.  SAR  data  and 
model  simulation  show  an  abrupt  increase  in  wavelengths 
between  some  of  the  solitary  wave  trains.  This  suggests  that 
the  model  formalism  does  contain  dynamics  similar  to  the 
ocean.  The  modeled  wavelengths  are  within  a  factor  of  2,  or 
less,  of  wavelengths  derived  from  SAR  data.  The  simulated 
phase  speeds  range  from  0.73  m/s  to  0.83  m/s.  The  phase 
speeds  estimated  from  the  measurements  range  from  0.8  m/s 
to  1.1  m/s. 

4.  Conclusion 

[53]  This  work  establishes  the  generation  and  propagation 
of  solitary  waves  in  the  Shandong  Peninsula  region  of  the 
Yellow  Sea.  It  is  shown  that  generation  and  propagation  of 
internal  solitary  waves  occurs,  over  small  topographic 
drops,  along  a  southeastern  track  off  the  Chinese  coast 
located  south  off  the  Shandong  peninsula.  SAR  imagery 
indicates  the  presence  of  internal  solitary  waves  along  the 
same  track  and  suggests  their  generation  in  the  shallower 
shelfbreak  region.  Model  results  confirm  the  generation  of 
internal  solitary  waves  in  the  same  off  shelf  area.  On-shelf 
SAR  imagery  and  model  results  show  the  presence  of 
internal  bores.  The  SAR  observations  are  used  for  guiding 
and  tuning  the  model  simulations  by  comparing  spectra  of 
observed  and  modeled  wavelengths.  The  tuned  model 
yielded  wavelengths  that  are  within  2  or  less  of  SAR  data. 
Modeled  wavelengths  ranged  from  335  m  to  3300  m  and 
amplitudes  up  to  9  m.  Modeled  phase  speeds  tended  to  be 
at  the  lower  limit  of  phase  speeds  deduced  from  SAR 
data  that  ranged  from  0.8  to  1.0  m/s.  Both  model 
simulations  and  SAR  observations  show  the  presence  of 
two  types  of  internal  solitary  wave  structures  configura¬ 
tions,  A/B.  This  suggests  that  the  model  formalism  does 
contain  dynamics  similar  to  the  ocean.  In  configuration  A, 
the  wavelengths  decrease  toward  the  back  of  the  packet. 
In  configuration  B,  they  increase  toward  the  back  of  the 
packet.  Configuration  B  tends  to  develop  later  in  time. 
Mixtures  of  both  configurations  can  occur  in  an  internal 
solitary  wave  train. 
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